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Abstract 

A mapping F2 population from the cross 'Piel de Sapo' x PI1 241 1 2 was selectively genotyped to study the genetic control of 
morphological fruit traits by QTL (Quantitative Trait Loci) analysis. Ten QTL were identified, five for FL (Fruit Length), two for 
FD (Fruit Diameter) and three for FS (Fruit Shape). At least one robust QTL per character was found, flqsS.l (LCD = 16.85, 
R^ = 34%), fdqsU.l (LCD = 3.47, R^=11%) and fsqsS.l (LQD= 14.85, R^ = 41%). flqs2.1 and fsqs2.1 cosegregate with gene a 
(andromonoecious), responsible for flower sex determination and with pleiotropic effects on FS. They display a positive 
additive effect (a) value, so the PI 1241 12 allele causes an increase in FL and FS, producing more elongated fruits. Conversely, 
the negative a value for flqsS.l and fsqsS.l indicates a decrease in FL and FS, what results in rounder fruits, even if PI1241 12 
produces very elongated melons. This is explained by a significant epistatic interaction between fsqsl.l and fsqsS.l, where 
the effects of the alleles at locus a are attenuated by the additive PI 1241 12 allele at fsqsS.l. Roundest fruits are produced by 
homozygous for PI124112 at fsqsS.l that do not carry any dominant A allele at locus a (PiPiaa). A significant interaction 
between fsqsS.l and fsqsll.l was also detected, with the alleles at fsqsll.l producing more elongated fruits. fsqsS.l seems 
to be allelic to QTL discovered in other populations where the exotic alleles produce elongated fruits. This model has been 
validated in assays with backcross lines along 3 years and ultimately obtaining a fsqsS.l -H\L (Near Isogenic Line) in 'Piel de 
Sapo' background which yields round melons. 
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Introduction 

Melon {Cucumis melo L., 2n = 2x = 24) is one of the most 
important horticultural crops, ranking in the ninth position of the 
world production of vegetables, with more than 27 million tons in 
2012 [1]. In the last two decades, melon production has 
experienced a 29% rise, what has been accompanied by an 
increase of 8% in the area harvested. It does not seem feasible to 
achieve the same gain in production in the next decades by 
increasing the land dedicated to its cultivation, so the efforts of 
plant breeders should be directed to improve traits having a great 
potential to increase yield and fruit quality. In this sense, fruit 
morphology can be considered a quality trait as it is one of the first 
features to be perceived. Although consumers associate particular 
sizes and forms to certain market classes (i.e. medium size and 
round cantaloupes or large and oval inodorus), a whole world of 
possibilities is being opened up before us with the development of 



products with original appearance, either by traditional breeding, 
like the ZESPRI GOLD kiwifruit cultivar [2], or by genetic 
engineering, like the purple tomato [3]. Modern market trends 
demand smaller melons, as the family size is progressively falling; 
and with rounder shapes, what minimizes costs and mechanical 
damages during storage and transport. 

Melon germplasm shows an impressive array of fruit sizes and 
morphologies [4-6] , ranging from a few tens of grams up to 10 Kg 
and from slightly flat to extremely elongated shapes. Most of this 
diversity is found in traditional Central Asian landraces and wild 
melons [7-8]. This gene pool is considered "exotic" when 
compared to the popular Occidental varieties (i. e. Galias, 
honeydews, western shippers, cantaloups, and inodorus like Tiel 
de Sapo', from here on referred to as 'PS') although it represents a 
valuable source of variability with a huge potential in breeding 
programs. In this sense, the accession PI124112 (C. ynelo ssp. 
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Figure 1. Fruit Shape (FS) polymorphism at each stage of the crossing program to obtain the fsqs8. 7-NIL. Fruits of (A) 'PS'; (B) PI 1241 12; 
(C) 2M158-3 (F,); (D) different F2 plants to show the transgressive segregation of FS; (E) 3M70-47 (Fj); (F) 5M113-1 x 'PS' (BC1); (G) 6M59-13 x 'PS' 
(BC2); (H) 7M36-1 x 'PS' (BC3); (I) 8M42-24 x 'PS' (BC4); (J) 9M7-15 : (BC4S1); (K) 10M2-30 x 'PS' (BC5); (L) 11M27 : (BC5S1), (L.I) 11M27-20: 
heterozygous at fsqsS.l and (L.2) 11M27-11: homozygous at fsqsS.l; (M) the fsqs8.1-N\l 12M57-3 OP (BC5S2). 
doi:10.1371/journal.pone.0104188.g001 



agrestis var. momordica) collected in 1937 in India is genetically 
very distant from commercial cultivars, as revealed by studies 
carried out with SSR (Simple Sequence Repeat) [5], [9], and 
AFPL (Amplified Fragment Length Polymorphism) and SNP 
(Single Nucleotide Polymorphism) [10] markers, and shows 
resistance to several fungal and viral diseases (reviewed in [7]). 

A few major genes with pleiotropic effects on melon fruit shape 
(FS) are currently known, like andromonoecious (a), that has been 
recently cloned [11] and is responsible for sex determination 
together with gynoecious (g), and like pentamerous ip), that controls 
carpel number. However, most of the phenotypic variation 
observed in FS seems to be caused by other Quantitative Trait 
Loci (QTL), as demonstrated by a number of QTL mapping 
studies tiiat have been undertaken to genetically dissect its 
molecular basis [12-18]. 

Genotyping technologies have experienced an authentic revo- 
lution in the last few years [19-20], what has led to a renewed 
interest in QTL mapping. The availability of affordable high- 
throughput genotyping platforms will enable researchers to 
annotate more precisely the genomic regions of interest and speed 
up drastically the formerly time-consuming process of obtaining 
valuable plant material. AU this will make possible to validate QTL 
and study putative interactions among them, at the same time that 



favorable alleles are transferred into elite varieties with a minimal 
drag of non-targeted fragments with potential detrimental effects. 
Consequendy, it is expected that in a near future, new QTL wiU 
be added to those already cloned in important crops, like the QTL 
for fruit sugar content Brix9-2-3 [21], for fruit size fw2. 2 [22], and 
those for fruit shape ovate [23], fasciated [24] and mn [25] in 
tomato; a QTL for grain shattering sh4 [26] and those for 
flowering date Hdl [27], Hd6 [28], HdJa [29], and EMI [30] in 
rice; and the QTL for plant architecture Thl [31] and flowering 
time Dwarfs [32] and Vgtl [33] in maize, among others (reviewed 
in [34—35]). In most cases, the strategy followed was the positional 
cloning, for what the availability of Near-Isogenic Lines (NILs) was 
essential, highlighting the importance of generating the appropri- 
ate genetic material to be successful. In this sense, the development 
of a collection of NILs in melon between the Korean cultivar 
'Songwhan Charmi' and TS' [36] has allowed to dissect in single 
Mendelian factors economically important traits as fruit morphol- 
ogy [14], [17], climacteric ripening [37] and virus resistance [38]. 

In the present work, we identify and localize several QTL 
associated with melon fruit morphological traits in a cross between 
the Spanish cultivar 'PS' and the Indian accession PI124112 and 
study the interaction between those with the strongest and most 
robust effects. Furthermore, we have validated the QTL for FS 
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fsqsS.l in a QTL-NIL developed in the genetic background of the 
cultivar 'PS', what produces totally round melons instead of the 
typical oval ones. 

Materials and Methods 

All experiments were conducted in compliance with the current 
Spanish laws. No specific permissions were required for the 
experiments carried out in these two locations, Cabrils and 
Zaragoza, as they are research institutes, Institut de Recerca i 
Tecnologia Agroalimentaria (IRTA) and Centre de Investigacion 
y Tecnologia Agroalimentaria de Aragon (CITA), respectively. 
Moreover, authors IE and JMA belong to IRTA and CITA, 
respectively. 

Plant material and phenotypic evaluation 

The F2 plants of the mapping population were derived from a 
single cross between the monoecious Indian landrace accession 
PI 1241 12 (momordica group) and the andromonoecious Spanish 
cultivar 'PS' (inodorus group) (Figure 1). PI124112 typically 
produces medium size, elongated and non-sweet fruits, whereas 
fruits coming from 'PS' are big, oval-shaped and very sweet [4—5]. 

A multi-location assay was conducted in summer of 2003 with a 
total of approximately 400 F2 plants. On one hand, 200 of them 
were randomly distributed in a greenhouse at IRTA in Cabrils 
(Barcelona, Spain), grown in drip-irrigated peat bags and self- 
pollinated by hand. The remaining 200 plants were randomized in 
an open field without irrigation at CITA, located in Zaragoza 
(Spain) and allowed to open-pollinate. The fertilization method 
became irrelevant as all the fruit except the seed content is derived 
from the maternal genitor. Ten plants of every genitor plus the Fj 
were included in both trials. 

Fruit length (FL) and maximum fruit diameter (FD) were 
measured on longitudinal sections of collected fruits and then, FS 
was calculated as the ratio of FL over FD. 

The absence of statistically significant genot)'pe-by-environment 
(GxE) interactions was tested by calculating the least square 
means of FS in a two-way ANO VA with the data coming from six 
to nine replicates ofFl, PI 124112 and 'PS' in Cabrils and 
Zaragoza as a necessary condition to merge the data coming from 
the two locations. 

Pair-wise Pearson's correlation coefficient was calculated 
between the three traits with the data coming from the whole Fg 
population and also with those obtained in the two location assays 
(Cabrils and Zaragoza) independentiy. 

Percent mid-parental heterosis (MPH) was calculated using the 
data derived from a total of 15 replicates of each of the genotypes, 
Fi and the both parents, as: 

MPH = {2FI -Fl -F2) / {Pi +P2) x 100 

where Fl, PI and P2 are the mean values of the Fi, 'PS' and 
PI 124 112, respectively, for the trait under study. 

One-way ANOVA and t-Student mean comparison tests were 
carried out with the data coming from the whole Fg population 
and also with those obtained in the two location assays 
independently. 

All statistical analyses were performed using JMP v5.1.2 
software for Windows (SAS Institute Inc., Cary, NC). 



DNA marker analysis in the mapping population 

A selective genotyping strategy [39], which exploits maximum 
phenotypic dilferences with a minimal labor and time investment, 
was approached to study the FS trait. The plants selected to be 
genotyped were those with values more extreme than the mean of 
the F2 population ±1SD from each location, resulting in 84 
individuals from both locations (File SI). 

DNA was isolated from young leaves according to [40], with 
slight modifications and genotyped with a total of 133 markers (66 
SSRs, 65 SNPs and 2 Cleaved Amplified Polymorphic Sites, 
CAPS) chosen to cover all the genome among those available in 
the literature (File S2). SSR markers were amplified following the 
protocol with M13-tailed forward primers plus dye labeled Ml 3 
primers described by [41] and subsequentiy analyzed on a 3130 
Genetic Analyzer (Applied Biosystems). Genotyping with SNP 
markers was carried out using the Sequenom MassArray iPLEX 
system (lUumina), followed by MALDI-TOF mass spectrometry at 
Servicio de Investigaciones Biomedicas, Unidad Central de 
Investigacion (University of Valencia, Valencia, Spain). Amplifi- 
cations of tiie CAPS markers, AluICAPS [1 1] and PSI_25-H03 
[42], were performed in 25-|j.l volume solutions containing 20 ng 
of genomic DNA, 7.5 mM Tris HCl (pH 9.0), 5 mM KCl, 2 mM 
(NH4)2SO,, 2.5 mM MgCl2, 0.2 mM of each dNTP, 0.2 jxM of 
reverse and forward primers and 2.5 U of DNA polymerase 
(Biotools, Madrid, Spain). PCRs were carried out on a TC-512 
thermocycler (Techne, UK) programmed with an initial denatur- 
ation step at 94°C for 5 min, followed by 35 cycles of 94°C for 
30 s, 62°C for 45 s and 72°C for 1 min, plus a final elongation 
step at 72°C for 5 min. A 10-|xl volume of the AluICAPS and 
PSI_25-H03 PCR products were digested with the restriction 
enzymes ^&I and fii/j 12361 (Fermentas, Fisher Scientific, Madrid, 
Spain), respectively, following the supplier's recommendations and 
separated by agarose gel-electrophoresis. 

Linkage map construction and QTL analysis 

The linkage map was constructed using MAPMAKER 3.0 [43] 
and the distances were calculated with the Kosambi function [44] . 

Markers were associated with the group command with LOD> 
3.0. When possible, they were put in order within eadi group 
using the order command with LOD>3.0, followed by a 
verification with the ripple command and the building of the 
framework map. If necessary, remaining markers in the group 
were included with the place command. If the previous approach 
did not render any map, the compare function was used to 
determine the best possible order for five starting markers, adding 
the remaining ones with the try command. In the case of groups 
with a low number of markers, the order and distances were 
estimated using the group command, followed by compare, what 
generated the framework map. Finally, maps were drawn with 
MapChart version 2.2 [45]. 

QTL analysis was performed for each location independentiy 
and for data from both locations merged together using Windows 
QTL Cartographer 2.5 [46] with the composite interval mapping 
(CIM, [47]) procedure. The LOD score threshold for a significant 
level p<0.05 was calculated by a permutation test with 1000 
resamplings, giving a value of LOD>2.5 for the whole population 
for all traits and largc'r numbers in the individual populations. 
Additionally, the threshold LOD > 2. 5 was also established for the 
single location analyses. Markers associated to the most robust 
QTL (6 out of the 10 QTL as several of them share the same 
marker as the strongest linked one, i.e.jlqsS.l and fsqsS.l) were 
genotyped in all F2 plants in order to obtain unbiased estimates of 
the gene action (d/[a]). The QTL were named according to the 
nomenclature system adopted by [48], where the first letters 
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designate the trait abbreviation, followed by a "Q," that stands for 
QTL, then a letter referred to the mapping experiment ("s", in this 
case), followed by a digit representing the hnkage group (LG) to 
which the QTL maps, and then followed by a dot and a final digit 
that distinguishes different QTL from the same experiment on the 
same LG. In those cases in which two overlapping QTL were 
identified in a limited region and for the same trait, they were 
effectively considered as only one QTL. 

Study of QTL interactions 

Epistatic interactions between the FS QTL detected were 
investigated by two-way ANOVA with theJMP v5. 1.2 software for 
Windows and using the genotypic data from the linked markers in 
all F2 plants (not only the subset employed for th(; selective 
genot\ping). These markers were AluICAPS in the case- o(fsqs2.1, 
CMAT141 for fsqsS.l, and ECM 67 (or fsqslZ.l. hitcraitions 
were considered to be statistically significant when j!)£0.05. Sum of 
squares values were used to calculate the percentage of the 
phenotypic variation explained by the interaction between two 
QTL (R2). 

QTL validation by developing a QTL-NIL 

The strategy followed to validate one of the FS QTL detected 
and to introgress it into the 'PS' elite cultivar consisted in selecting 
a single F2 plant from the trial located in Cabrils, that carried the 
PI124112 alleles at the QTL region and also showed round fruit, 
and backcross it to the parental TS' (Figure SI). Five rounds of 
backcrossing guided by a doubled selection, phenotypic (favoring 
round-shaped fruits) and assisted by markers linked to the target 
QTL (to monitor the introgressions harboring PI 1 24 1 1 2 alleles), 
were carried out. Simultaneously, background markers covering 
all the genome were also employed to gradually increase the 
percentage of the recurrent parent genome by selecting against the 
donor (PI 124 112) genome. An advanced backcross selection, 
10M2-30, was genotyped with 548 SNPs using the Golden-Gate 
assay (lUumina, San Diego, CA) described in [10] in order to verify 
introgression size and integrity, and the composition of the 
genomic background. After two cycles of self-fertilization, the 
definitive introgression line in the BC5S2 generation was obtained 
(Figure IM). 

Phenotypical evaluations of the fruits were conducted in the 
greenhouse at several stages of the backcross scheme. In 2008, 42 
fruits collected from ,37 plants belonging to three full-siblings BG4 
famiUes (8M42-37, 41 and 42) carrying the PI1241 12 allele for the 
FS QTL in LG VIII (homozygous and heterozygous for two linked 
markers placed 23.3 cM apart in the genetic map) were evaluated 
in Cabrils. In 2009, BC4S1 plants (Figure SI) were genotyped with 
the markers linked to the QTL and 30 of their fruits were 
phenotyped in a greenhouse at Polytechnic University of Valencia. 
The cosegregation of markers and FS phenotypes was assessed by 
ANOVA. FinaUy, in 2012, melon fruits from the BC5S2 family 
12M57 and the parental cultivar TS' were analyzed in Paiporta 
(Valencia). 12M57, derived from a homozygous plant for the 
closest markers to the FS QTL fsqsS.l (11M27-11, Figure SI), 
contained the introgression for the PI124112 allele at those same 
markers linked to the QTL. Means of the fruits produced by 10 
12M,57 plants were compared with means of the fruits collected 
from 15 TS' plants by a t-test. 

Results 

Phenotypic variation in morphological traits 

Cultivar TS' produces oval-shaped fruits (FS=1.4±0.1; 
Table 1 and Figure lA), while PI124112 displays elongated fruits 



(FS= 1.7±0.2; Table 1 and Figure IB). The Fi exhibited higher 
FS values (FS = 2.0±0.2, Table 1 and Figure IC), in accordance 
with best-parent heterosis for this character previously observed in 
melon [16], particularly in hybrids between exotic accessions and 
TS' {inodorus group), caused mainly by an increase in FL [5]. No 
significant GxE interaction was observed for FS, therefore the 
parental mean values were calculated with the data coming from 
both locations. In Cabrils, the MPH values were 77%, 10% and 
60% for FL, FD and FS, respectively; in the assay located in 
Zaragoza, 120%, 34% and 68%, respectively; and for the whole 
set of data, 100%, 24'X) and 64"'(i. rc'spcrti\-cly. In all cases, 
heterosis was confirmed and, though there are slight difierences in 
the MPH values for each character among assays, the ranking 
remains unaltered. 

A continuous distribution of values was observed in the F2 
(Figure 2), as expected in a quantitative trait like FS, with values 
ranging from 0.86 (oblate fruits) to 3.23 (extremely elongated 
fruits), making clear a transgressive segregation for this character 
in this population (Figure ID and Figure 2C). Transgressive 
segregation was also observed for FL. In both cases, it was 
observed in both directions in the whole F2 population (Table 1 
and Figure 2), although the phenotypes where more extreme for 
the higher values, that means fruits with a more elongated shape. 
In all cases, individuals with fruits exceeding the mid-parent value 
±3.84SD were found. The trait showing a broader range of values 
was FL (from 7.1 to 38 cm), displaying extreme values corre- 
sponding to —6.56 and 4-13.52 times the SD of the mid-parental 
FL. Conversely, FD showed the lowest variation (from 7.30 to 
17.70 cm), what represents the mid-parent FD value — 3.84SD 
and +5.55SD, respectively. In between these two cases, the FS 
values ranged from — 4.31SD and +11.23SD the mid-parent 
values, respectively. 

FL and FS showed a highly significant (p<0.001) and strong 
correlation with values for the Pearson's correlation coefficient of 
0.84 for the whole population, and 0.86 and 0.83 for the 
subpopulations of Cabrils and Zaragoza, respectively. In contrast, 
FD and FS were not significantiy correlated in any case. 

Linkage map construction 

One hundred and twenty-two out of the 1 33 total markers (63 
SNPs, 57 SSRs and 2 CAPS) used to gen()t)'pc the F2 mapping 
population could be incorporated to a geni^tic map consisting in 18 
LGs that span 1437.3 cM (Figure 3), what renders a density of 1 
marker every 11.8 cM. By comparing with the melon consensus 
map [48] , each LG could be identified, finding that five of them 
were partial LGs of the integrated map. The map shows an 
average marker gap of 13.8 cM. 

QTL associated with fruit morphology traits 

Given that no GxE interaction was found among locations, 
data from both locations were merged for subsequent QTL 
analysis. A total of ten QTL for fruit morphology-related traits 
were positioned along the melon genetic map (Table 2 and 
Figure 3) with the selective genotyping strategy, five underlying FL 
(in LGs II, Illb, Via, VIII and Xb), two responsible for FD (in LGs 
Ilia and XII), and tiiree controlfing FS (in LGs II, VIII and XII). 
LOD score values ranged from 2.80 to 16.85, in the case of 
fdqsja.l and flqsS.l, respectively. The QTL with the highest 
effect for each of the three characters studied were flqsS.l 
(LOD = 16.85, = 34%), fdqsl2.1 (LOD = 3.47, R^ = 1 1 %) and 
fsqsS.l (LOD = 14.85, R^ = 41 %). The remaining QTL accounted 
for 4 to 17% of the phenot)'pic variation. LOD peaks of the two 
major QTL located in LG ll,flqs2.1 andfsqs2.1, and in LG VIII, 
flqsS.l and fsqsS.l, overlapped completely (Figure S2). 
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Table 1. 


Mean and Standard Deviation (SD) values of the traits for the parents, the F, 


and the whole Fj. 








FL 




FD 




FS 






Mean 


SD 


Mean 


SD 


Mean 


SD 


'PS' 


19.5 


1.3 


14.0 


1.1 


1.4 


0.1 


Pn24112 


14.8 


1.8 


9.2 


1.1 


1.6 


0.2 


F, 


25.1 


3.5 


12.7 


1.8 


2.0 


0.2 


F2 


19.5 


5.7 


11.7 


2.0 


1.7 


0.5 



FL: Fruit Length; FD: Fruit Diameter; FS: Fruit Shape. 
doi:l 0.1 371 /journal.pone.01 041 88.t001 



The interval defined by the markers flanking the QTL ranged 
from 12.9 (flqs3h.l) to 69.1 (flqs6a.l) cM, being less than 26 cM 
for five out of the ten QTL (Table 2). 

The QTL on LG VIII (flqsS.l and fsqsS.l) were consistentiy 
detected in the three analyses (that is, with the merged data and 
with those commg from individual locations, File S3). The QTL 
fsqs2.1 was also detected in the three analyses with a 
LOD = 2. 5. The rests of the QTL were identified in two of the 
analyses, except tor flqslOb.l . In all cases, the additive efTects sign 
remained invariable across the difTerent environments. 

Since fsqs2.1, fsqsS.l and fsqsl2.1 were the QTL most 
consistently detected among analyses, all plants were genotyped 
with markers AluICAPS, CMAT141, and ECM67, located in the 
respective QTL region, in order to obtain unbiased estimates of 
their effects (Table 2). Same direction effects of the exotic allele 
(PI 124 112) for FL and FS QTL with the same chromosomal 
location were observed. In this sense, flqs 2.1 and fiqs2.1 displayed 
a positive value for the additive effect (a), 3.06 and 0.27, 
respectively, what means that the PI124112 allele causes an 
increase in FL and FS, producing more elongated fruits. 
Conversely, the negative value of a for flqsS.l and fsqsS.l 
(—4.95 and —0.50, respectively) is indicative of a decrease in both 
characters, FL and FS, what results in rounder fruits, even if 
PI124112 itself produces very elongated melons. In all the 
remaining QTL, except for fdqsl2.1, detected with the data 
coming from Cabrils and with the whole set {a = — 1.03 and a = 
— 1.10, respectively), the PI124112 allele contributed to increase 



the value of the particular morphological trait (Table 2 and File 
S3). So, the PI124122 allele atfdqsl2.1 will render elongated fruits 
by decreasing their diameters, whereas in the case oiflqs2.1 and 
fsqs2.1, a similar phenotype will be obtained by increasing the 
length of the fruit. 

Gene action values ranged from —0.35 and —0.32 (flqsS.l and 
fsqsS.l, respectively) to 0.75 flqs2.1 andfsqs2.1; Table 2). In most 
cases, the alleles at the QTL loci showed an additive effect, except 
for flqs2.1 and fsqs2.1, where the inheritance mode is near 
complete dominance. In general, the gene action values obtained 
for the same QTL using the subsets of data coming from both 
locations (File S3) agree with the ones calculated when all the data 
were pooled together, with slight differences that can be attributed 
to the smaller sample size of the first analysis. 

The additive effect at fsqsS.l is supported by the phenotypic 
data, as the heterozygous l^YL-fsqsS.l displays a FS value (i.e. 
1.24) intermediate between the two homozygotes, the homozygous 
line for the PI 1241 12 allele (i.e. 1.06) and the 'PS' parent (average 
FS= 1.40). 

QTL X QTL interaction 

In order to detect the putative effects of other genomic regions 
on each of the FS QTL, all the possible digenic interactions 
between them were studied although only those that were 
statistically significant when /;<0.05 are shown (Table 3 and 
Figure 4). 
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Figure 2. Frequency distribution of the traits across tlie F2 population derived from the cross 'PS' x PI1 241 1 2. (A) FL (Fruit Length), (B) 
FD (Fruit Diameter) and (C) FS (Fruit Shape). Both parents and F, values are marked. In the upper box, 25*, 50* and 75* quartiles are displayed; the 
sample mean and the 95% confidence interval are represented by a diamond, and the outlier values as dots. The bracket along the edge of the box 
stands for the part of the graph in which 50% of the observations are gathered together. Mean and standard deviation estimates for a normal 
distribution are also shown. 
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Figure 3. Location of 10 QTL responsible for morphological fruit traits mapped in a F2 population derived from the cross 'PS' x 
PI1241 12. FL: Fruit length, open bars; FD: Fruit Diameter, hatched bars; and FS: Fruit Shape, solid bars. Length of QTL bars corresponds to the two- 
LOD support intervals (from peak marker) based on the results obtained by CIM using Windows QTL Cartographer 2.5. Numbers on the left side 
correspond to the distance in cM (according to [44]) from the top of each chromosome. 
doi:1 0.1 371/journal.pone.01 041 88.g003 



fsqs2.1 X fsqsS.I 

A significant interaction was identified between fsqs2.1 and 
fsqsS.l (p = 0.013, Figure 4A, Table 3), that explains a 4% of the 
genetic variance of FS. 

flqs2.1 has been previously identified as gene a [12], [49] and 
the marker linked to it (AlulCAPS) is actually the causal mutation 



of andromonoecy [11]. For that reason and from now on, the 
alleles at this locus will be referred to as A and a. 

It is well known that the dominant allele A causes elongated 
shapes in some genetic backgrounds [12], as becomes clear in the 
cases of homozygous for PS allele at locus fsqsS.l (PsPs). For each 
of the three possible genotypes at locus a (AA, Aa, aa), the lowest 
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Figure 4. Digenic interactions studied by two-way ANOVA between tlie FS QTL using the mariners significantly linked to them. (A) 

fsqsS.l and fsqslA (fsqslA x fsqsS.l); (B) fsqsTlA (fsqsS.l x f5qsl2.1). Alleles at locus flqslA are named A and a since this QTL has been previously 
identified as gene a [12], [49]. PsPs: homozygous for the allele PS (solid line); PsPi: heterozygous (dashed line); PiPi: homozygous for the allele 
PI124112 (dotted line). 
doi:1 0.1 371 /journal.pone.01 041 88.g004 



values of FS index (rounder fruits) were always observed in the 
homozygous for the PI124112 allele at locus fsqsS.l (PiPi) 
(Figure 4A). So, when the PI124112 allele is present at fsqsS.l 
and the dominant A allele is absent at locus a (PsPiaa and PiPiaa), 
there is a considerable decrease of FS. 

In the homozygous PiPi, the differences among the three 
dififerent genotypes at locus a are the slightest. It seems that 
homozygous PiPi at fsqsS.l diminish or, in some way, attenuates 
the effects of the dominant A allele of gene a on fruit elongation. 

In the case of heterozygous genotypes at fsqsS.l (PsPi), some 
differences can be observed depending on the alleles present at 
locus a (A_ are expected to produce moderately elongated fruits 
while aa will rendered more rounded ones). So, we can only see 
marked differences in the three possible genotypes at gene a when 
the PI 1241 12 allele is not present at QTh fsqsS.l (no masking 
effect) or it is in heterozygosis, PsPs or PsPi (moderate attenuation 
of the phenotype caused by A), with FS showing nearly parallel 
increments when these two scenarios are compared. In both cases, 
the elongated phenotype caused by the aforementioned dominant 
A allele becomes clear, as homozygous AA produce the most 
elongated fruits, followed by the heterozygous (Aa). 

fsqsS.l X fsqs12.1 

A significant interaction was identified between fsqsS.l and 
fsqsl2.1 (p = 0.014, Figure 4B, table 3), explaining 7% of the 
genetic variance of FS. 

Table 3. Significant digenic interactions between all the FS QTL. 



Similarly to what can be observed in the interaction y5^,s2. 7 x 
fsqsS.l , the lowest values of FS (rounder fruits) are obtained when 
the PI124112 aUele is in homozygosis at fsqsS.l QTL (PiPi) 
(Figirre 4B). In such a case, the phenotypic variations due to the 
different allelic composition at fsqsl2.1 locus are practically 
imperceptible. Those differences in the FS value become more 
evident with the presence of the PS allele M fsqsS.l, specifically 
when it is in homozygosis (PsPs). Interestingly, the highest value of 
FS (the most elongated fruits) are reached for the PS homozygous 
at fsqsS.l locus that are homozygous for the PI1241 12 at fsqsl2.1 
(PsPsPiPi), what is totally in line with the positive additive value of 
fsqsl2.1 QTL (the Pi allele causes an increase in the FS index). 

Here, unlike what happens in the previous interaction, within 
the same genotype for QTh f sqsS. 1 , the increase in FS due to the 
PI 1241 12 aUele at fsqsl2.1 becomes clear when the PI124112 
allele is absent at locus fsqsS . 1 . 

QTL validation in the fsqsS.l -NIL 

The identification and mapping of the fsqsS.l QTL in the Fj 
population has allowed us to use marker-assisted selection to 
precisely introgress the PI 1 24 1 1 2 allele responsible for the round 
shape into the elite cultivar TS' as the result of a backcross 
program (Figure SI). To ensure the integrity and maintenance of 
the QTL effects, two rounds of QTL verification were carried out 
in intermediate generations. 



QTL pair 


Source 


DP 


Sum of Squares 


F Ratio 


Prob> F 


fsqs2.1 -fsqsS.l 


fsqs2.1 (AIulCAPS) 


2 


7.83 


35.72 


<0.0001 




fsqsS.l (CMAT141) 


2 


26.79 


122.21 


<0.0001 




fsqs2.1 (AIulCAPS) x fsqsSl (CMAT141) 


4 


1.42 


3.24 


0.0132" 


fsqsS.l-fsqsU.l 


fsqsSl (CMAT141) 


2 


23.02 


82.76 


<0.0001 




fsqsl2.1 (ECM67) 


2 


1.59 


5.71 


0.0038 




fsqsSl (CMAT141) x fsqsl2.1 (ECM67) 


4 


1.78 


3.19 


0.0142" 



The interactions were explored by two-way ANOVA using the genotypes for the closest markers in the whole F2 population derived from the cross 'PS' x PI124112. 
": Interactions were considered to be statistically significant when p<0.05. 
doi:l 0.1 371/journal.pone.Ol 041 88.t003 
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In 2008, the effect of the fsqsS.l was monitored with markers 
CMAT14I and PSI_25-H03 in BC4 plants (Table 4, 8M-codes in 
Figure SI). Plants carrying the PI1241 12 allele in heterozygosis at 
fsqs8. 1 produced significantly rounder melons than plants carrying 
the PS allele in homozygosis (F = 21.10 and F = 25.94 for 
CMAT141 and PS_25-H03, respectively, p<0.001). 

In 2009, the three genotypes axfsqsS.l (monitored with marker 
CMAT141 and PSI_25-H03) segregated in a BC4S1 population 
(9M7, Figure SI), re-confu~ming the cosegregation of PI124112 
alleles at fsqsS.l with the round fruit phenotype (Table 4; 
F= 13.75 and p<0.001 for GMAT141; F = 5.71 and p<0.01 for 
PSI_25-H03). 

At the next stage, whole-genome genotyping was performed in a 
selected BC5 plant (10M2-30) in order to check the integrity of the 
introgressed region and the composition of the genetic back- 
ground. The target introgression harbored 8.5 Mb in the LG VIII, 
and there were only seven additional introgressions ranging from 
250 kb to 86.8 Mb. 

Finally, the definitive introgression line was frxed in a BC5S2 
family (12M57, Figure SI). The famUy 12M57 showed highly 
significant shorter and rounder fruits, that is, lower values of FL 
and FS, respectively (Table 4; F= 10.27 and F= 105.71, for FL 
and FS, respectively, and p<0.001). Additionally, significant 
effects on FD were also observed (F= 10.90, p<0.001. Table 4). 

Discussion 

The Fg population between TS' and PI124112 showed an 
impressive variability for fruit morphology. Transgressive segre- 
gation in both directions was found for FL and FS, indicating that 
at least one of the parents contributed alleles with opposite effects 
and making possible to obtain genotypes that could render round 
fruits. In fact, the existence of an allele inducing round-shaped 
melons was identified in the genitor PI1241 12 and verified with an 
introgression line in the TS' background. The effects of this allele 
seem to be masked by epistatic interactions with other gene(s) in 
the PI124112 background, and those effects were unveiled in the 
TS' background. A similar discovery has been reported recentiy 
with the obtaining of a climacteric line derived from the cross 
between two non-climateric melon accessions, TS' and 'Shong- 
wan Charmi' [37]. This is a further indication of the richness of 
the melon gene pool that remains stiU largely unexplored. 

The segregating population coming from parents genetically 
distant and with morphological differences used in this work has 
revealed itself as suitable for the genetic dissection of morpholog- 
ical traits and mapping QTL. Furthermore, the most robust FS 
(and the highly correlated FL) QTL detected (flqs2.1, fsqs2.1, 
flqsS.l and fsqsS.l) were consistent across the locations tested, in 
agreement with other results published [16], in which QTL 
affecting FS displayed a high heritability, what explained their 



reproducibility in assays carried out in different locations and 
along several years. 

The availability of a melon consensus map including the 
position of QTL controlling a panoply of important traits [48] 
made easier the comparison of the QTL identified here with those 
previously described. In this sense, most QTL controlling the same 
characters studied in this work have been detected previously and 
mapped to similar genomic positions (Table 2). In fact, homolo- 
gous QTL obtained in different populations can be assigned to all 
the QTL described here, except in the case of flqslOh.l (only 
detected with the full set of data) and fdqsBa.l (not found in 
CabrHs). Only a QTL associated to FL was recently reported in 
the LG II of the related cucurbit Zucchini [Cucurbita pepo L.) in a 
region syntenic to the LG III of melon at which fdqs3a. 1 maps 
[50]. However, this QTL could not be confirmed in further 
backcross generations. 

The QTL in LG II (flqs2.1 and/or fsqs2.1) was previously 
identified as a pleiotropic effect of the gene a [12], [49], whose 
alleles are responsible for sex determination in melon flowers [51] 
as it has been reported in other plant species like Vitis ssp. [52]. 
The dominant inheritance mode observed for the alleles Mflqs2.1 
and fsqs2.1 loci also supports it. The cloning of the gene has 
vmveiled that the causal mutation of the andromonoecious melon 
types is in the active site of 1 -aminocyclopropane- 1 -carboxylic acid 
synthase [11]. 

The map position oi fsqsS.l (and flqsS.l, as they overlap) 
happened to be similar to those reported for other FS QTL 
(Table 2; reviewed in [48], [bmfsqcS.S {8.1 in [14]) and fsqnS.l 
(fshS.l in [54]). Those QTL were detected in a NIL population 
derived from the cross TS' (C. melo var. inodurus) x PI161375 
'Shongwan Gharmi' (C. melo var. chinensis), and a Recombinant 
Inbred Line (RIL) population coming from the cross PI4 14723 (C. 
melo var. momordica) x 'Dulce' (C. melo var. reticulatus), 
respectively. So, fsqsS.l seems to be aUehc to those major QTL 
underlying fruit morphology in a broad array of melon 
germplasm. The PI161375 and PI414723 alleles produced 
elongated fruits. On the contrary, PI 1 24 1 1 2 allele (in a certain 
allelic combination in trans with gene a variants) induced round 
fruits, what may be showing the existence of an allelic series at this 
QTL when different germplasm sources are investigated. Three 
more QTL have been reported in the surrounding region (fsqaS. 1 , 
fsqnS.l and fsqcS.l, compiled in [48]) and two additional ones 
further away in the same chromosome (fsqcS.l and fsqpS.2, 
reviewed in [48]). These results suggest that an important part of 
the loci controlling FS in melon could be concentrated in the LG 
VIII. It is also possible that some of them are the same QTL that 
the one described here (harboring common and/ or new alleles). 
To discern between these two cases, higher resolution in the 
mapping and QTL studies is required. Similar QTL clusters have 



Table 4. F-values obtained from ANOVA to validate the effects of QTL fsqsS.l on FS. 







2008 




2009 




2012 


CIV1AT141 


PSI_2S-H03 


CIV1AT141 


PSI_25-H03 


CIV1AT141 


FL 


5.58' 


8.76" 


10.10"' 


4.82* 


10.27"' 


FD 


3.46™ 


2.49"' 


0.53™ 


0.38™ 


10.90" 


FS 


21.10™ 


25.94"" 


1 3.75"' 


5.7 1" 


105.71'" 



Data coming from the assays carried out along 3 years were used. Plants were genotyped with linked markers (CMAT141 and/or PSI_25-H03}. 
FL: Fruit Length; FD: Fruit Diameter; FS: Fruit Shape. 
'"P<0.001;"P<0.01; "P<0.05; ™P>0.05. 
doi:l 0.1 371/journal.pone.Ol 041 88.t004 
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Eilso been observed for other traits in tomato [55] and maize [56] - 
[57]. Furthermore, all the QTL for FL, which is strongly 
correlated to FS (the authors didn't consider the FS ratio as a 
trait), identified by [58] in cucumber (another member of the 
genus Cucumis), were located in LG 4, what is syntenic to LG VIII 
of C. melo together with LG 6 of cucumber [59-60]. In agreement 
with this, most of the FS QTL found by [61] were located in 
cucumber in those two same chromosomes, 4 and 6. In Zucchini, 
however, the reported major QTL involved in fruit shape was 
mapped in LG VI, which is syntenic to LG IV of melons and Chr3 
of cucumber [50]. 

If most of the fruit morphological variability present in 
cultivated melons is attributable to an array of alleles at a limited 
number of loci, as it seems happens in tomato [62]), then their 
unambiguous characterization provides scientists and breeders 
with an useful tool to engineer it. 

The FS QTL fsqsl2.1 identified in this work maps in the 
neighborhood of the chromosomal segment harboring the gene p, 
which controls carpel number. Pleiotropic effects of p on fruit 
shape has been reported previously [12-13], similarly to what 
happens with the gene a. Also, fruit shape variation has been 
related to differences in the number of carpels in C. pepo. A QTL 
associated to this variation has been mapped in LG V in a region 
of C.pepo genome that is syntenic to that one in melon LG XII in 
which the p gene is located [50] . 

In a series of NILs derived from 'PS' (produces melons with 
three carpels) and the accession 'Shongwan Charmi' (renders fruits 
with five carpels), fruits with five carpels were rounder than those 
with three [14], [16-17]. However, a high number of carpels is an 
undesirable feature in commercial types as it causes the fruit to 
have a bigger empty cavity. The carpel number did not segregate 
in the population studied in this work, therefore the variation in 
the number of carpels is not a consequence oifsqsl2.1 effects. On 
the other hand, it is also possible that at least part of the fruit 
morphology variation associated to p in the other populations 
could be caused by another linked FS QTL. The fsqcl2.1 QTL 
was mapped within a 5 cM interval by [17], further genetic 
resolution would be necessary to discard the presence of another 
linked QTL. 

When QTL responsible for FS identified in tomato, eggplant 
and pepper have been compared, an interesting scenario is 
revealed, as they could be orthologous loci controlling the same 
character in all these species. In line with this, the same connection 
perhaps can be established with melon. To do it, the recendy 
obtained melon genome sequence [60] would be of great value. 

Positive values of the additive effect of QTL in LG II, flqs2. 1 
and fsqs2.1 (Table 2), indicate that alleles from PI 124 11 2 act in 
the direction of increasing the magnitude of both FL and FS. The 
same can be argued for fiqsl2.1 QTL regarding FS. Surprisingly, 
it is quite the opposite for the QTL in LG VIII, flq.sS.l and 
fsqsS.l, which display negative additive effects (Table 2). So, for 
this QTL, the exotic PI124112 alleles contribute to decrease the 
values of FL and FS, rendering rounder fruits even if PI1241 12 
itself produces very elongated melons [5] . Furthermore, it is well 
known that the dominant alleles of gene a (A), responsible for the 
monoecious sex type (A_), also causes elongated melon fruits in 
most genetic backgrounds. However, andromonoecious cultivars 
(aa) that produce elongated fruits are also common (i.e. TS'). 
Evidence for very significant interactions between these two QTL, 
fsqs2.1 and fsqsS.l (Table 3 and Figure 4A), may shed Ught on 
this. So, if no interaction would exist, round fruits would be only 
produced by homozygous for the exotic allele at fsqsS.l that do 
not carry any dominant A allele at locus a (PiPiaa). It seems that 
there is a requisite that must be satisfied in order to obtain a round 



fruit and that is to have two copies of the PI1241 12 allele at the 
fsqsS.l locus. Therefore, only the homozygous (PiPi) can 
potentially render rounded melons. Furthermore, for these 

homozygous PiPi, the presence of the A allele at gene a has 
practically no effect on the elongation of the fruit and it is then 
when the epistasis between these two loci becomes obvious. Only 
in this scenario, the real effect of the PI1241 12 alleles at fsqsS.l is 
clear. Among the heterozygous at fsqsS.l, only those not having 
any A allele at locus a (PsPiaa) will be able to produce moderately 
elongated fruits. The rest (PsPiAA and PsPiAa), will give very 
elongated fruits. According to this, the homozygous PsPsAA 
generates the most elongated fruits, followed by PsPsAa. 

The additive effect of alleles at fsqsS.l has been confirmed by 
the phenotype of the heterozygous NlLrfsqsS.l , with a FS value 
intermediate between the two homozygotes. This agrees with the 
results observed by [16] in a population derived from the cross TS' 
X 'Shongwan Charmi', in which the QTL alleles inducing fruit 
enlargement were dominant, whereas those inducing rounder fruit 
were additive or recessive and only minor epistatic interactions 
were found. 

So, the action of the QTL fsqsS. 1 attenuates the effects of the 
alleles at locus a. Therefore, the dififerences due to the genotype at 
the latter are more evident when the PI124112 allele is absc-nt at 
fsqsS.l (PsPs, in our case), as it has been previously described in 
cases of QTL epistatic interactions in other fruits, like tomato [63]. 

Adding complexity to this scenario, the other interaction 
detected between two FS QTh, fsqsS.l and fsqs 12.1 , .shows a 
similar pattern though, in this case, the exotic allele M fsqsl2.1 
causes the enlargement of the fruit. This means that the roundest 
fruits will be produced by individuals that are homozygous for the 
PI124112 allele fsqsS.l and for the 'PS' aUele &t fsqs 12. 2 
(PiPiPsPs). On the other hand, the most elongated fruits will be 
obtained from the other possible double homozygous (PsPsPiPi). In 
between, the heterozygous wiU show a wide range of FS values, 
with the double heterozygous exhibiting an intermediate pheno- 
type. 

To conclude, the effect of the PI124112 allele tor fsqsS.l 
remains unnoticed in the PI124112 accession. Now, this becomes 
understandable because alleles for FS QTL with opposite effects 
have been detected in the accession PI 124 112. To unveil the 
interactions reported here has been possible thank to the isolation 
of the effect of fsqsS.l from the interference of other FS QTL in 
the fsqs8.1-mh. 

The literature is riddle with examples of epistatic interactions 
among QTL affecting the same character, like ljfs2 (probably 
corresponding to ovate), ljfs7 (likely being sun and having 
pleiotropic effects on pericarp elongation), ljfs3 (it maps close to 
fs3.2) and jfs 11 (in the same region that fsll.l), all of them 
affecting a morphological trait like eccentricity in tomato [64]. A 
strategy successfully employed to uncover QTL hidden by epistatic 
interactions with other genes has been to cross a certain QTL-IL 
(introgression line) with different genotypes, so new phenotypes 
arise in certain genetic backgrounds [63]. 

Conclusions 

The additive inheritance mode of fsqsS. 1 , with the exotic allele 
coming from PI124112 being the one causing the round shape, 
make the fsqsS.l -NIL a more suitable material to modify the 
morphology of commercial varieties by obtaining hybrids. 
Actually, monoic cultivars normally produce elongated fruits and 
for that reason, they are not commonly used to generate new 
commercial hybrids as the FL and FS values increase as a result of 
heterosis [5], [16]. However, monoecy brings some advantages 
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from a breeding perspective since it eliminates the need for hand 
emasculation during the crossing process, and specifically to obtain 
hybrids, that usually produce fruits with improved quality traits. 
Therefore, unveiling and understanding epistatic relationships 
between the QTL affecting FS and gene a, like the one discovered 
with fsqsS.l, could make possible the development of hybrids that 
bring together both desirable characteristics, that is, potential new 
monoecious varieties producing round-shaped fruits. 

Supporting Information 

Figure SI Pedigree showing the time scale and the crosses 
carried out to generate the mapping F2 population and to 
eventually obtain the/i^iS.i-NIL. 
(TIF) 

Figure S2 Locations of QTL for FL (Fruit Length) and FS (Fruit 
Shape), /Z552.i and fsqs2.1, respectively, in chromosome II; and of 
jlqsS.l .und fsqsS.l on chromosome VIII, showing their overlap- 
ping. The threshold LOD (2.5) was calculated for a significant 
level p = 0.05 by a permutation test (n=: 1000). Peak marker: in 
bold; flanking markers: underhned. 
(TIF) 
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